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ABSTRACT 

This  study  deals  with  the  unsteady  free  convective  transient  flow  of  an  incompressible  dissipative  viscous  fluid 
past  an  infinite  vertical  plate  in  presence  of  viscous  dissipative  heat,  under  the  influence  of  a  uniform  transverse  magnetic 
field.  The  velocity  profiles,  temperature  profiles,  skin  friction  and  rate  of  heat  transfer  have  been  discussed  with  the  help  of 
graphs.  The  problem  is  governed  by  a  coupled  non-linear  system  of  partial  differential  equation.  The  exact  solutions  are 
not  possible  hence  explicit  finite  deference  method  is  employed. 

KEYWORDS:  Non-Linear  System,  Velocity  Profiles,  Field  of  Engineering 

INTRODUCTION 

The  problem  of  free  convection  and  mass  transfer  flow  of  viscous  fluid  under  the  influence  of  magnetic  field  has 
attracted  the  interest  of  many  authors  in  view  of  its  application  to  geo-physics,  astrophysics  and  engineering.  Boundary 
layer  control  in  the  field  of  aerodynamics  and  so  on.  The  researches  on  flow  of  viscous  fluids  are  being  used  in  almost  all 
field  of  Engineering  including  Bio-medical  Science,  meteorology,  physical  chemistry,  plasma  physics,  Monsoon  dynamics 
and  allied  field.  The  unsteady  conducting  free  convection  flow  of  dissipative  fluids  past  an  infinite  plate  have  received 
much  attention  because  of  non  linearity  of  the  governing  equations,  without  taking  into  account  viscous  dissipative  heat 
and  magnato-hydrodynamic.  This  problem  was  discussed  by  Siegal  (1958)  by  integral  method.  Some  of  the  literature 
surveys  and  review  of  pertinent  work  in  this  field  are  documented  by  Goldostein  and  Eckert  (1960),  Sparrow  and  Gregg 
(1960),  Chung  and  Anderson  (1961),  Gebhart  (1961),  Menold  and  Yang  (1962),  Schetz  and  Eichhorn  (1962)  and 
Goldestien  RJ.  and  Briggs  (1964),  In  all  these  problems  the  effects  of  viscous  dissipative  heat  and  Magnetohydrodynamic 
was  assumed  to  be  neglected  but  Gebhart  (1962)  has  shown  that  when  the  temperature  difference  is  small  or  in  high 
Prandtl  number  fluid  or  when  the  gravitational  field  is  of  high  intensity,  viscous  dissipative  heat  should  be  taken  into 
account  in  steady  free  convection  flow  past  a  semi-infinite  vertical  plate.  The  following  assumption  Soundalgekar,  Bhatt 
and  Mohiuddin  (1979)  discussed  the  effects  of  free  convection  current  on  the  flow  past  an  impulsively  started  infinite 
isothermal  vertical  plate.  The  free  convection  and  mass  transfer  effects  on  the  flow  past  an  infinite  moving  vertical  porous 
plate  with  constant  suction  and  heat  sources  when  free  stream  velocity  is  an  oscillatory  function  of  time  was  presented  by 
Raptis  (1982). 

The  combined  buoyancy  effect  of  thermal  and  mass  diffusion  of  magneto  hydro  dynamic  natural  convection  flow 
was  presented  by  Agrawal,  Ram  and  Singh  (1980).  The  free  convection  heat  transfer  between  two  long  vertical  plates 
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moving  in  opposite  direction  was  presented  by  Kajravelu  (1978).  Recently-  Soundalgekar,  Lahurikar  and  Pohanerkar  have 
studied  free  convection  flow  of  an  incompressible  viscous  dissipative  fluid;  this  problem  is  governed  by  coupled  non  linear 
system  of  partial  differential  equation  solved  by  finite  difference  technique.  Present  study  is  unsteady  free  convection  flow 
of  a  dissipative  viscous  fluid  past  an  infinite  vertical  plate,  on  taking  into  account,  the  viscous  dissipative  heat  under  the 
influence  of  uniform  transverse  magnetic  field.  This  problem  is  governed  now  by  non-linear  coupled  differential  equations 
whose  exact  solution  is  not  possible.  So  we  employ  explicit  finite  difference  method  for  its  solution.  In  present  study  we 
consider  the  unsteady  transient  free  convection  flow  of  an  incompressible  dissipative  viscous  fluid  past  an  infinite  vertical 
plate  on  taking  into  account  viscous  dissipative  heat  under  the  influence  of  a  uniform  transverse  magnetic  field.  The 
velocity,  temperature  coefficient  of  skin  friction  rate  of  heat  transfer  and  the  effect  of  magnetic  field  parameter  M,  the 
Prandtl  number  Pr,  and  the  Eckert  number  E  have  been  discussed  with  the  help  of  graph. 

LIST  OF  SYMBOLS  USED 


x\y’ 

Cartesian  coordinate  system 

u 

Velocity  of  fluid 

g 

-A 

Acceleration  due  to  gravity 

K 

-> 

Thermal  conductivity  of  fluid 

p 

-A 

Density  of  fluid 

V 

Kinematics  Viscosity 

P 

-A 

Coefficient  of  volume  expansion 

Ten 

-4 

Temperature  of  plate 

Too 

-A 

Temperature  of  fluid  for  away  from  plate 

P 

—> 

Viscosity  of  fluid 

a 

-A 

Electrical  conductivity  of  fluid 

pe 

Magnetic  permeability 

Cp 

-A 

Specific  heat  at  constant  pressure 

U„ 

Reference  velocity 

TR 

— > 

Reference  time 

L 

Reference  length 

0 

-A 

Dimensional  temperature 

Pr 

Prandtl  number  parameter 

E 

-A 

Eckert  number  parameter 

M 

Magnetic  field  parameter 

q 

-A 

Heat  per  unit  area  per  unit  time 

11 

Dimensionless  parameter 

HO 

-A 

Magnetic  field  strength 

FORMULATION  OF  THE  PROBLEM 


Consider  the  unsteady  hydrodynamic  free  convection  flow  of  an  incompressible,  viscous  liquid  past  an  infinite 
vertical  plate.  At  initially  the  temperature  of  the  plate  and  the  fluid  are  assumed  to  be  the  same.  The  time  t'  >  0,  the  plate 
temperature  is  raised  to  tco  which  is  then  maintained  constant  under  these  conditions.  The  flow-variable  are  functions  of  f 
and  y'  alone.  The  x1  -axis  is  taken  along  the  plate  in  the  vertically  upward  direction  and  the  y'-axis  is  taken  normal  to  it.  A 
uniform  magnetic  field  of  intensity  HO  is  applied  in  the  y-direction.  Therefore  the  velocity  and  the  magnetic  field  are  given 


by  *7  (lEO)  anc|  H  (0, H0)  rp^g  ||ujt|  being  slightly  conducting  the  magnetic  Reynolds  number  is  much  less  than 
unity  and  hence  the  induced  magnetic  field  can  be  neglected  in  compression  with  the  applied  magnetic  field  in  the  absence 
of  any  input  electric  field  as  Sparrow  E.M.  and  Cess  R.D.  (1972).  Then  the  governing  equation  of  flow  are- 
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du'  d2u 

~W~vd^ 


P 


„  dT'  „  d2T'  ( du  A 
pC—  =  K—  +  M 


dt' 


dy 


dy. 


boundary  condition  are  - 


t'<  0, 

u  =  0, 

rp\ _  p\ 

for 

all  y' 

t'>  0, 

u'=  0, 

rT~1 1 _  rT~'  1 

1  ~  1  CO 

at 

y' 

u'=  0, 

^  —  *  CO 

at 

y1— >  °o 

(1) 


(2) 


(3) 


we  introduce  the  following  non-dimensional  variables 


-2 
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_  'T’t  yi 

~  1  1  G)' 

U0  =  (vg/3AT)U3 , 

L  = 

fgjBAT) 

2 

3 

9 

T  = 

1R 

( gjBAT ^ 

3 

K  v  ) 

l  v  J 

t' 

y-£ 

u  =  - 

u' 

0_  (T'-T'J 

tr 

L 

U  o 

pCp 

u 2 

E  -  U°  . 

M 

 °p:ulrR 

K 


CpAT 


Using  the  above  variables  of  equation  (4)  equations  (1)  and  (2)  reduced  to- 

du  d2U  n  .. 
dt 


P'f  =  0+Pr£ 


f duV 


dy 


the  boundary  condition  (3  )reduce  to  - 

t<  0,  u  =  0,  9  =  0  for  all  y' 

t>  0,  u  =  0,  9  =  1  at  y  =  0  > 

u  =  0,  9  — >  0  at  y  — >  <*> 


SOLUTION  OF  THE  PROBLEM 


(4) 


(5) 


(6) 


A  new  dimensionless  variable 


9  = 


y 

i+y 


In  view  of  (8)  the  equation  (5)  to  (7)  as 


(8) 
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du 

dt 


(l-?7)4 


-2(1  -r/f 


d u 

—  +  6 -Mu 
dr/ 


(9) 


Pr— =  (l-77)4^-2(l-77)3  — +  (l-/7)4Pr  E 


dt 


dr/ 


dr/ 


du 

dr/ 


t<  0,  u  =  0,  0  =  0 
t>  0,  u  =  0,  0  =  1 

u  —  0,  0  — >  0 


for  all  r/ 
at  r/  =  0  > 
at  r/  — > 


(10) 


(ll) 


Here  equation  (9)  and  (10)  are  non-linear  coupled  partial  differential  equations  and  are  to  be  solved  by  using  the 
initial  and  boundary  condition  (11).  However,  exact  solutions  are  not  possible  for  this  set  of  equations  and  hence  we  solve 
these  equations  by  finite-difference  method.  The  equivalent  finite  difference  scheme  of  equation  for  (9),  (10)  and  (11)  are 
as  follows: 


U‘j+1  UU  =d  -%j) 

At 


4  j  Ui+lj  2Uij  +  U,_hJ  1  .3  J  i+l.j  ^  ij 


(A  r/f 


2(1  -r/uY 


A  r/ 


■7+1  wiJ  -n_„  x4  J  20 \j  +  Qt_x  J 


Pr^+I  °i  J  =  (1 


At 


(A  r/Y 


■2(l-r/uf\di+Li 


A 1/ 


+(i-^)4p4^H 


(12) 


(13) 


Here  index  i  refers  to  17  and  j  to  time.  The  mesh  system  is  divided  by  taking  Ar|  =0.1  from  the  initial  condition  in 

(11)  as 

w(  0,0)  =  0,  0(0,0)  =  1,  u(i,  0)  =  0,  0(z,  0)  =  0 

for  all  i  except  i  =  0 

w(0,  j)  =  0,  0(0,  j)  =  l  for  all  i 

u(\,  j)  =  0,  0(1,  j)  =  0  for  all  i  (14) 

First  the  velocity  at  the  time  step  viz.  u(i,j+  1)  (i=I,  10)  is  computed  from  (12)  in  term  of  velocity  and  temperatures 
at  points  on  the  earlier  time  step.  Then  0(i,j+l)  is  computed  from  (13).  The  procedure  is  repeated  until  t=l,  (i.e.  j=800). 
During  computation  At  was  chosen  as  0.00125.  These  computations  were  carried  out  for  Pr(=0.71,  7)  and  E(=0,  0.1,  0.3, 
0.5)  to  judge  the  accuracy  of  the  convergence  and  stability  of  finite-difference  scheme,  the  same  program  was  run  the 
smaller  values  of  At,  i.e.  At  =  0.0009,  0.001  and  no  significant  change  was  observed.  Hence  we  conclude  that  the  finite- 
difference  scheme  is  stable  and  convergent. 


NAAS  Rating:  2.73-  Articles  can  be  sent  to  editor@impactjournah.us 


Unsteady  Free  Convection  Transient  Flow  o  f  an  Incompressible  Dissipative  Viscous  Fluid 


13 


SKIN  FRICTION 

Given  non-dimensional  form  as 


T  =  ■ 


du 

dr] 


T  =  - 


PU20 


^ =°  here 

Numerical  values  of  r  are  calculated  by  applying  the  Newton's  interpolation  formula  for  four  points. 

HEAT  TRANSFER 

The  non-dimensional  form,  the  rate  of  heat  transfer  can  be  shown  to  be  given  by 

'  do 

1  j(7=0 


dr/ 


We  have  computed  q  by  following  the  above  procedure.  These  values  of  q  are  listed  in  table  1. 

RESULT  AND  DISCUSSION 


The  velocity  profit  (u)  versus  (q)  for  different  values  of  magnetic  field  parameter  M  and  time  t  respectively  for 
cases  of  the  Prandtl  number  Pr  (=0.71  and  7.00).  Then  we  see  that  the  velocity  profile  (u)  increases  with  increase  in  time  t, 
where  as  velocity  (u)  decreases  with  the  increase  in  M,  in  all  the  cases  Pr  (=0.71  and  7.00).  Hence  we  conclude  that 
velocity  (u)  decrease  with  increase  in  Prandtl  number  also  we  observed  that  velocity  (u)  increases  first  near  the  plate  and 
then  trend  gets  reversed  as  q  increase  which  is  shown  in  Fig  (1)  and  (2).  The  Temperature  distribution  0  is  drawn  against 
17  for  different  value  of  Pr  and  1.  We  see  that  0  decreases  with  increase  in  Pr  {(3),  (4),  (5)}.  Further,  we  observe  that  0 
increase  in  t  {(I),  (2),  (3)},  here  also  we  conclude  that  0  decrease  with  increase  in  17  which  is  shown  in  Fig  (3).  The  skin- 
friction  is  shown  against  t  for  different  value  of  M  in  the  cases  of  Pr  (=0.71  and  7.00)  respectively.  Here  we  observed  the  T 
decreases  with  increase  in  M  or  Pr  where  as  it  increases  with  the  increase  in  t  which  is  shown  in  Fig  (4).  We  observe  that  q 
increases  with  the  increases  in  M  where  as  it  decreases  with  increase  in  t.  Also  we  observe  that  u  is  drawn  against  q  for 
different  value  of  Eckert  number  E  in  this  case  velocity  (u)  increases  with  increase  E  for  all  cases  of  Prandtl  number. 
Further  we  observe  that  q  increases  with  the  increase  in  Pr  which  is  shown  in  table  1 . 
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Figure  1:  Velocity  Profile 


Figure  2:  Velocity  Profile 
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Figure  3:  Temperature  Profile 


Figure  4:  x  Against  t  for  Different  M 


Table  1:  Value  of  q  when  E=0.3 


Pr 

t/M 

1.0 

5.0 

7.0 

0.71 

0.1 

1.46691025310025 

1.466952321643101 

1.46696932150880 

0.2 

1.0263125091201 

1.02634225069120 

1.02640350921620 

0.3 

0.60213401961801 

0.60314231629002 

0.68201509213321 

7.0 

0.1 

4.72351692031120 

4.72352931462902 

4.72382931624123 

0.2 

3.02134296132013 

3.02136923497213 

3.02213416923812 

0.3 

2.01132945612312 

2.01142936916654 

2.01182345618213 

100.0 

0.1 

15.93215628932126 

15.932172931461123 

15.93221345216113 

0.2 

14.90235436891234 

14.902376542131234 

14.90247532162214 

0.3 

13.30506921354278 

13.30508963125981 

13.30507921462293 
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